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ABSTRACT 


Photon absorption coefficients and mean opacities have been calculated for 
hydrogen, beryllium, carbon, nitrogen, aluminum, and silicon over a tempera- 
ture range from 1.5 to 34 ev and a aensity range from about 107} g/cm? down- 
ward. Contributions to the absorption coefficient from free-free (inverse- 
bremsstrahlung), bound-free (photoelectric), and bound-bound (line-absorption) 
processes are included, as is Compton scattering. Certain thermodynamic 
properties are also given. An improved recipe for pressure ionization is de- 
rived which is approximately valid at nondegenerate densities for any ratio of 
Debye length to ion- sphere radius. Line absorption is evaluated using recent 
results from pressure-broadening theory and a representation of line series 
which is computationally as simple as the statistical method. The results show 
that lines increase the Rosseland mean opacity by a factor which can be nearly 
ten and which is insensitive to moderate changes in line widths. The code 
employed generates ionic energy levels internally by isoelectronic interpolation, 
and is immediately applicable to any mixture of elements in which no ion has 
more than 14 bound electrons. 

The results of the calculations of thermodynamic properties and mean 
opacities are given in the tables in Appendix A, and the graphs of the mono- 
chromatic absorption coefficients are given in Appendix B, which comprises 
Vols. II and III. 
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I. BASIC DEFINITIONS AND RELATIONS 


When photons travel in an absorbing medium, the probability that a 
photon goes a small distance, dx, without interacting with the medium is 
1 - dx/A, where A is called the absorption mean free path. The reciprocal 
of A is called the linear absorption coefficient, cross section per unit 
volume, or inverse mean free path, and is conventionally denoted by yp. 
The ratio p/p, where p is the mass density of the medium, is the mass 
absorption coefficient or cross section per unit mass, k. It is the sum of 
the absorption cross sections of a representative ensemble of the atoms, 
ions, and molecules in the medium, divided by the sum of their masses. 
Generally, A, #, and « depend strongly on photon frequency and thus carry 
the subscript v. Certain weighted averages of se with respect to v are 


called opacities. In particular, the '"Rosseland mean opacity" is definedas 


ie dB 
dv 
0 aT 


*R co dB (1) 
v 1 q 
[ dT ~hv/kT . 
0 -e % 


cp 2 —__,——__ (2) 


where 


hv/kT aged 


B= B (T) = Gir hey te 1) 


Vv 


is the specific intensity (energy per unit time, area, solid angle, and 


frequency interval) at frequency v of black-body radiation at tempera- 


1 


2 


ture T, and where 


00 4 
B= | B dv= ot 
Jo." id 


The opacities which have just been defined are relevant to radiative 
transfer problems only when the matter is in local thermodynamic equilib- 
rium (LTE) at temperature T, and in the case of the Planck mean this 
condition is not always easy to establish. When it holds, the spontaneous 
emission coefficient (energy radiated per unit time, mass, solid angle, and 


frequency interval) is 


; -hv/kT 

j,F Be (l-e ) (3) 
whether the radiation present is in equilibrium with the matter or not. If 
the radiating region is now assumed to be optically thin (small compared 
with the photon mean free path, A) at all frequencies, every photon emitted 
will escape, there is no induced emission, and the total rate of energy loss 


per unit time and mass is just 47x_B; but this overestimates the energy- 


loss rate if any photons are ee i.e., if the radiation present 
approaches black-body intensity at any frequencies. Accordingly, the 
utility of the Planck mean is restricted to situations where particle collisions 
are frequent enough to maintain LTE without any assistance from the radi- 
ation field. 

The Rosseland mean is relevant only to an optically thick region in 
which a small temperature gradient exists--more precisely, in which 
A, |VTI << T and AY is much less than the dimensions of the region, at all 
frequencies. Then the radiation field at any point looks nearly like the 
black- body radiation appropriate to the local temperature (except for a 
small asymmetry in angle, which drives the radiative flux), and LTE is 
maintained automatically, whether particle collisions or radiative processes 


dominate the transition rates between states. Under these circumstances, 


the radiative flux (net energy per unit time and area crossing a surface 


normal to the temperature gradient) is 


A 
a R 
F = -41 =~ VB, (4) 
where AR = 1/pk,. Derivations of the above relations are given, for 


(1) 


example, by Chandrasekhar. 
There is a large class of radiative transfer problems involving 
regions which are optically thick at some frequencies and thin at others, 
so that even if LTE prevails, the Rosseland and Planck mean opacities are 
not applicable. An example is the atmosphere of almost any star. In such 
problems it is generally necessary to proceed iteratively, treating each 
frequency separately and using the full frequency-dependent absorption 
coefficient. 
A process which enters the radiative transfer problem in a somewhat 
different way is scattering, in which a photon changes direction but not 


(2)(3) 


frequency. It is readily shown that the effect of scattering on the 


radiative flux, when the Rosseland mean is applicable, is correctly included 
if in Eq. (1) (1 - eee, is replaced by (1- aad | + o where 
ie (1/2) ies (1 - to (t) dt and a(t) is the cross section per unit mass 
for scattering through an angle arc cost. When the Planck mean is appli- 
cable, scattering should not be included, since it is not an energy-loss 
mechanism. 

In much of the present work, we use dimensionless variables, which 
we define here. In place of the photon frequency v, we use the reduced 


frequency variable 


_ hy 
eins Ss ues 
in terms of which 
BY dv B. du 15 ae 
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The combination Ky (1 - ) which appears in the Rosseland and Planck 


means is actually the absorption coefficient corrected for induced emission 


and is often called Ki This gives the compact forms 


foe) CO 
+ = | Why) du ; Kp = | k' b(u) du . 
R O°. SG oe, . 


Since part of the absorption coefficient varies inversely with the cube of the 


frequency, we also use the quantity 
3 
D(u) = i 
(u) =u K 


and the Stromgren function 


so that 
foo) foo 
1 1 dS 15 -u 
7 meres (ene du . K = — ! D(u)e du . 
Kp i D(u) du Pp nt 0 


In the specification of the thermodynamic state of plasmas, it is most 
convenient to use as independent variables the thermal energy, kT (in 
electron volts when not otherwise stated), and the dimensionless degeneracy 


parameter 


a 2(2nmkT)>! 2 
ies i cae eae 
h N 
e 
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(where No is the number density of free electrons and m is the electron 
mass), which is, roughly, the phase space (in units of h?) available per 
free electron, or the factor by which the existing electron density would 
have to be increased isothermally before strong degeneracy effects would 
set in. In all our applications, Tis muth larger than unity and the energy 


distribution of free electrons is nearly Maxwellian. 


ar cai 


II]. OCCUPATION NUMBERS AND THERMODYNAMIC PROPERTIES 


The calculation of the populations of the various stages of ionization 
and excitation in a plasma is a prerequisite to finding its radiative proper- 
ties. The thermodynamic properties, which are readily obtained fromm the 
relative populations, are also intrinsically useful. Our treatment of this 
part of the problem follows standard lines, except for the pressure- 
ionization recipe used. We start, for a pure element, with a set of energy 
levels, Es statistical weights, gi and net ionic charges, fi. The Boltzmann 
and Saha equations give the thermal-equilibrium unnormalized population, 


P.. of state i: 


BE, 4-1 
i 
P, =e, explaf, - i - is) ; (5) 
where @ = 4n land Je is the pressure lowering of the ionization potential 
(in units of kT) of the ion species of charge f (thus, Jo refers to the first 
ionization potential). Jy is evaluated in Section III. The normalized popu- 


lations, or occupation numbers, Pp,» are given by 


Pp. = , (6) 


and all of the thermodynamic properties are averages over the P;- Thus, 


we have the mean and mean-square charges, 
<— 2 2 
o>) = oS Pf; and (2 » a P,f; ‘ (7) 
i i 


the total internal energy per nucleus, 


> = 27 {B +e (f, +1) - a) (8) 
and so on. Since we are using kT and @ (or I) as independent variables, the 
electron density, N. and electron pressure, P. = N. kT, are known at the 
start, but the mass density, p, and total pressure, P, require a knowledge 
of (2): 

Das ie 1 
leery aa Poti h) Po (9) 
where A is the atomic mass number and M, is the unit atomic mass. 

The convergence of the partition-function sum 2), P. is secured by 
setting P. = 0 for all states whose loosest electron is bound (in the isolated 
ion) by less than the perturbation, JkT, of the ionization potential. As will 
be seen later, the existence of such a cutoff is necessary only for the 
thermodynamic part of the calculation; the contribution of the high bound 
states to the sum giving the absorption coefficient, once (a ‘has been 
determined, converges anyway. 

For preliminary estimates, the concept of the ''dominant potential" 


(4)(5) 


is useful. It is defined, in our notation, as akT, and ‘na thermal 
ensemble of ions the most abundant are those having ionization potentials 
close to the dominant potential. This permits (2) , etc., to be estimated 
immediately. Stated differently, the one-electron levels with binding energies 
larger than the dominant potential are those which, on the average, are 
occupied nearly to their full capacity. 

In our results we have not included the contribution of radiation to the 
pressure and internal energy. This contribution is simply additive; for 
equilibrium (black-body) radiation the energy density is 

4nB : 


© a “kT ergs 
ere ee Geren —) =e (10) 


and the radiation pressure is one-third of this, i.e., 


4 
_ 4"B _ kT dynes 
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Il]. PRESSURE IONIZATION 


In a material of finite density the potential distribution in and near a 
given ion is influenced not only by its own bound electrons, but also by free 
electrons, by neighboring ions, and (slightly) by neutral atoms. These 
perturbers produce effects of two kinds: their time-averaged effect is to 
alter the set of energy levels available to the ion in question, and their 
time-dependent effect is to broaden these levels both by shifting them back 
and forth adiabatically and by inducing transitions between them. The first 
effect, besides providing a natural cutoff to the bound-state partition func- 
tion, effectively lowers all the ionization potentials and shifts the equilibrium 
occupation numbers in the direction of increased ionization. This effect, 
which is most important at high densities, is usually called pressure ioni- 
zation. The second effect, pressure broadening, is discussed in Section VI. 

Pressure ionization has often been treated by assigning to each ion 
a sphere occupying the ion's share of the total volume and containing the 
ion and enough free electrons to make the sphere's net charge zero. In 
the simplest form, the free-electron density is assumed uniform. Inside 
the sphere the potential due to other spheres is neglected. Thus, close 
encounters between ions are disregarded, and the picture is basically that 
of a crystal lattice, with strongly correlated ion positions. This picture 
may be expected to hold at high densities and low temperatures. At the 
opposite extreme is the case of nearly random ion positions; the first-order 
deviation from randomness leads to the Debye-Huckel potential, which is 
valid in the low-density, high-temperature limit. In this section we develop 
a theory of pressure ionization which yields the ion-sphere and Debye-Huckel 
results as limiting cases. 

Consider a nucleus, Z, fixed in a sea of electrons and point ions at 


a temperature T. We wish to find the time-averaged potential distribution 


ll 


12 

around Z, and we assume for the present that only this average is felt by 
any one particle; that the potential distribution is spherically symmetric 
around Z; that it is determined by the time-averaged, spherically symmetric 
charge density; and that it changes only slightly within one particle wave- 
length. We describe the electrons by Fermi statistics and the ions by 
Maxwell statistics. These assumptions are those of the finite-temperature 
Thomas-Fermi (TF) atomic model, extended to include the neighboring 

ions. The electrostatic potential, $(r), is then determined by the Poisson 


equation 


az 


ae | Na 
—V o> = Fe (ro) = “ane 2 ain - Rey , 


where n, is the electron density, n. is the number density of ions of charge 


z., and 
i 


n,(r) = n, (co) exp (2, <2 ) : 


00 ,/2 ae 
Fin | ea o(00) = 0; ro~> Zeasr~+0. 
Ose 44 


The degeneracy parameter, @, is determined by 


3/2 
n (co) - £(émmkT) 2 F(-a) . 


h? No 
In our applications @ is always positive and appreciably larger than unity, 
i.e., the free electrons are nondegenerate and 


= 
Na 


Close to the nucleus, where ed > @kT, there is a region of degenerate 


F(-a) fe 7% <<1. 


bound electrons which is, in fact, the ion core. The asymptotic ion densi- 


ties n. (co) satisfy the condition of electrical neutrality, > : zn, (co) = n (00). 
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Under the substitutions 


dwe® ae 2 
= kT 2) (2, + 2,)n, (00) 


ms 
= 


(note that D is the Debye length including ions and electrons), the above 


equations take the nondimensional form 


& qi [EY -a)  @ oper) 
2 ’ 


ae Lea 


x | 


maees Sp a aotee an average, weighted with n. (co), over ion species, and 
are DILOy The second term on the aes side may be simplified with 
small loss of accuracy by noting that it is comparable with the first only 


for fairly small y; in this region, 


cepa Gla ie pian eee}, 


so that to this approximation the assortment of ion species present can be 
* 

replaced by a single fictitious species with charge z , which is never less 

than unity, even when <z) +0. Our basic equation for the potential distri- 


bution is then 


2 
ae (xy) = == Beenn - expl-2"y] (12) 
Zz 


with boundary conditions y(oo) = 0 and xy > Ze“ /(DkT) = Ky asx->0., 
Before proceeding to the solution of Eq. (12), let us consider what 
information we hope to extract from it. We will use the TF model solely to 
obtain the perturbing potential produced by free electrons and neighboring 
ions, and not to obtain all of the ionic properties. Since the density distri- 
bution of these perturbers depends somewhat on the bound-electron distri- 
bution, we must assume a bound-electron distribution in order to find the 
perturber distribution and its resulting potential. Having done so by means 


of the TF model, which approximates the average bound-electron distribution, 


14 


we may use this same perturbing potential in conjunction with better repre- 
sentations of the bound electrons, such as the experimental levels of 
isolated ions. 

The perturbing potential, oo satisfies a Poisson equation in which 
only the free-electron and ion densities appear, and in dimensionless form 


this equation is 


2 x 
ld _-_!1 /Fly-qy)_.-zy 
ae od Bea [ F(-a) | © I (13) 


where v = eo,/kT and Fly - a, y) = (ty s(etta-y +1) dt. (Note that 
F(y - a, y)/Fly - @) is the ratio of free-electron density to total electron 
density in the TF model.) Unlike y(x), which contains the nuclear potential, 
v(x) is finite at the origin, and.v'(0) = 0. Taking v(oo) = 0, v(x) will be 


negative; we let -v(0) equal J. Then 


° se 1 (* 2 
v(x) = -J + J S(t)t dt - — I S(t)t” dt, (14) 
0 490 


where S(x) is the right-hand side of Eq. (13), containing the solution y(x) 
6 
of Eq. (12). Hence, J =, . S(x)x dx. We can put Eq. (14) ina more 


convenient form, 


ee hoes 2 Ta) aes 2 i 
v(x) = - = i S(t)t dt + I S(t)t dt - S(t)t dt , (15) 
0 x 


x 


in which the first term, dominant at large x, may be identified as the 
asymptotic Coulomb potential of the total excess of free electrons over 
neighboring ions. It will be seen shortly that the total potential y(x) vanishes 
exponentially at large x. Now the potential due to the nucleus and bound 
electrons alone, y(x) - v(x), has an asymptotic Coulomb form dependent 

only on the net ionic charge; this permits evaluation of the first term of 

Eq. (15) without a knowledge of y(x) all the way to the origin, provided that 


the ionic charge is known. 


15 


In the solution of the foregoing equations, the following expansions 
are useful; they are valid when @ is somewhat larger than unity, so that 
e” >> 1 and the free electrons are nondegenerate (only this case is treated 
in what follows): 


ite ee pa <a) 


Ely - a) _ 
F(-a) a2 Oe oinat? scat pees. 
Pi Set. 2 f+ eV - ext vy) = (y >> 1) 
vo prey, 


ly <<). 


From these and the form of Eq. (12),we can distinguish four regions of 
progressively smaller x (larger y). 

Region A. The limiting form of Eq. (12) for large x, where y << 1/2", 
is 


a’ 


eo ay oe ee 
: xy) =y> 


& 


with the solution y = (C/x) e *, which is the form of the Debye-Huckel 
potential. In the region where this holds, the charge densities of ions and 
electrons nearly cancel and most of the electrons are free. 

Region B. Proceeding inward, if z* >> 1, there is a region where 
1/2* << y << land the density of the ions and that of the bound electrons are 
small compared with that of the free electrons, which is approximately 


constant and equal to its asymptotic value. In this region, Eq. (12) becomes 
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approximately 


Ts Gr ei 
x dx” as + 1 
with the solution 
2 
pS 4B —>—_, 
x 6(z +1) 


in which A and Bare arbitrary constants. This recalls the ion-sphere 
model with uniform free electrons. 

Region C. In the next region, 1 << y <<a, the bound electrons out- 
number the free, but the occupation of the bound states is well below their 
capacity and the local kinetic-energy distribution is still approximately 
Maxwellian. The free-electron density is now larger than its asymptotic 


value. Here, 


2 
= Se (xy) # — ey . 
Ase zt+i 


Region D. Finally, in the ion core, y >> @ and most of the electrons 
are in fully occupied states. As the nucleus is approached, the bound- 
electron density diverges as sae and the free-electron density as gr are 
(fictions peculiar to the TF mode}); the total potential, y(x), approaches 
K)/x + const. + Gt *y, and the perturbing potential, v(x), approaches 
-J + ox! 25, 

The procedure for obtaining v(x) should now be clear. There are 
basically three parameters in the problem: ny a, and Ko: More conven- 
ient than the last of these, for purposes of inward integration, is C, the 
normalization constant of the large-x solution. We can choose 2 a, and 
C, integrate Eq. (12) inward, and obtain v(x) from Eq. (15). The net ionic 


charge, z, then is given by 


eee a 2 
Bs de ore DET 2 K= (i S(t)t dt, 
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and we recall that 


co 
-v(0)= J = { S(t)t dt. 
0 


We have carried out a number of numerical integrations for various values 
of z a(22), and C. The results may be summarized as follows: 

1. Kand J are virtually independent of @ for fixed z* and C, even 
though Ky isnot. That is, as one might expect, the ionic charge 
rather than the nuclear charge determines most of the free- 
electron distribution. 

2. v(x) is nearly equal to its value at the origin, -J, through all of 
region D (the core) and much of region C. Thus, most of the 
bound states of the perturbed ion, including all the fully occupied 
ones in the core, are simply shifted relative to the continuum by 
the same amount; their wave functions are like those of the iso- 
lated ion and their ionization potentials are lowered by JkT. This 
result contrasts with that of Keller and Meyerott; (> 7) what they 
actually found was the outer screening by all the electrons--and 
the bound contribution is already present in the isolated ion. 

A very simple approximate solution for Eqs. (12) and (15) is readily 

obtained by assuming that only regions A and B need be considered, i.e., 


that all the bound electrons are at the origin and the free-electron density 


is approximately uniform. Then Eq. (12) is just 


2 y x 
ld ~(2*41 
= , (xy) = = {1 -e (z Yj ; 
dae zt+il 


For large z* we approximate the right side by 


y when (z* + ljyy <1 (region A) , 


* . 
w + > ; 
co hen (z ljy > 1 (region B) 


wettest bald 
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and letting x) denote the transition point, we require y and dy/dx to be 
continuous at xy At small x, y ~ (K/x) - J +--+ , where to this approxi- 
mation the first term is the Coulomb potential of the ion and the next is the 
screening potential evaluated at the origin, i.e., the reduction in ionization 
potential (in units of kT) produced by the free electrons and ions. Using 


the solutions in regions A and B, which are, respectively, 


2 e* (x > x_) 
x 1 
a K ae 
—- (x <x,), 
x 6(z* + 1) 1 
and the continuity conditions 
ee 
-x 
Ce t= K-JIx = 
1 6(z* + 1) 
x 
-x 
Ce } =J- —_—— ' 
2(z" + 1) 


we get 


a Ge 1 2 
eos F)-she [(x, + 1) -1j, 


2 
x x 
1 1 l 3 
K = L+x, + —) = ———  [(x, + ])) -1J2I5, 
is  ( 1 >) 3(z* +1) 1 


or, eliminating x 


]’ 


22° + 1)d +1 = [3(2* + 1)K+ ‘Odes (16) 


which for small (z* + 1)K becomes J = K, and for large (z™ + 1)K becomes 


J = (1/2(z* + 1)] [3(2* + KP! 3. Defining the ion-sphere radius, a, for 


an ion of charge z by the condition 4na?/3 = z/n (00) and recalling the 


19 


expression for the Debye length (D), 


1 ieee 


p* kT 


(z* + 1) n (00) ; 


we see that the parameter 3(z* + 1)K is just 


3(2* +1) ee . 32 7 eS 2 ion sphere volume 
DkT 4nD-n (co) D Debye sphere volume © 
e 


When this parameter is small, 


when it is large, 


ee 
2 akT 

These values of J are, respectively, the results of the Debye-Huckel and 
ion-sphere models. Very roughly, one may say that the model giving the 
smaller depression of the ionization potential is closer to the truth. 

The scaling in z* + 1 which appears above suggests plotting the results 
of the numerical integrations in the form (z* + 1)J versus (z* + 1)K, with 
z* asa parameter. When this is done, it is found that the curves for all 
values of z* give J-values which are (1. 06 #0. 08) times the J given by 
Eq. (16), provided that (z* + 1)K < 3 (as is the case for all of the values of 
a, T, and Z we are concerned with in this report). 

In our calculations of the equilibrium occupation numbers, we proceed 
as follows: for a given temperature and @, an initial set of occupation 
numbers is obtained from the Boltzmann and Saha equations, using ionization 
potentials lowered by an amount JkT calculated with z* = 1, and including 
all the ionic states bound by more than JkT. Then z* and the Debye length, 


D, are found, as is K for each z; this fixes J for each z. In each stage of 


ionization we then discard all bound states in which the loosest electron is 
bound by less than JkT in the isolated ion, and we lower the ionization 
potential by JkT; the new set of states and ionization potentials, inserted 
into the Boltzmann and Saha equations, gives new occupation numbers. The 
iteration converges rapidly and one cycle is usually enough. In the above, 
we take z as the ionic charge with the electron in question removed; for 
example, the first ionization potential is lowered by JkT calculated for z= 1. 
Numerical experiments have been performed which evaluate the sensi- 
tivity of the thermodynamic properties and of the opacity to J. The results 
show that changing all J's by 50% changes thermodynamic properties by 
~1% and opacities by~ 4%. For this reason, in all of our work we have 


simply used 1. 06 times the J given by Eq. (16). 


IV. CONTINUOUS ABSORPTION COEFFICIENTS 


The total photon absorption coefficient (not including scattering) at 
any frequency of interest in our temperature range has, in general, contri- 
butions from three types of elementary transitions involving electrons: 

(1) photoionization (bound-free), (2) inverse bremsstrahlung (free-free), 
and (3) line absorption (bound-bound). The cross sections for the first two 
types of transitions are said to comprise the ''continuous" part of the 
absorption coefficient; hydrogenlike expressions for them can be found in 
any text on stellar atmospheres or interiors (e.g., see Ghandeasekhar!” 


Aller, ae) Unséld, (9) and Schwareschiia”)), Straightforward but approximate 
generalizations to nonhydrogenic atoms and ions have been employed in all 
calculations of plasma opacities known to us, and ours will be no exception. 
Physically, it is assumed that all radial integrals which enter transition 
matrix elements can be adequately evaluated by using hydrogenic wave 
functions scaled to some effective charge that is the same in the initial and 
the final states of a given transi’ion. No matter how this effective charge 
is chosen, this assumption does violence to the frequency distribution of 
oscillator strength from certain initial states (an extreme example is the 
ground state of any neutral alkali); we feel that this is by far the most 
serious limitation on the over-all accuracy of our calculations. (Some 
compensation of errors occurs, of course, since the hydrogenic and actual 
oscillator strengths are subject to the same sum rule.) Calculations such 


(11) (12) who do not make 


as those of Varsavsky and of Burgess and Seaton, 
the equal-effective-charge assumption, show deviations from the scaled 
hydrogenic results for individual transitions which are large compared with 
the additional error made by setting the hydrogenic Gaunt factors equal to 
unity, which we have done throughout. 


Consider the photoelectric ejection of an electron of type j from an 
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atom or ion in state i, leaving a residual ion in state k. We choose the 
effective charge ms for this transition as determined by the ionization 
potential of electron j, evaluated without pressure-ionization corrections 
(this is because outer screening causes Z, to be an underestimate of the 
charge which determines the initial wave function, and pressure ionization 


effectively increases outer screening). This is the same choice used by 
(13) (14) 


Armstrong, Holland, and Meyerott and by Bernstein and Dyson. Thus, 
z 
= = 7 7 
7) Ry I. EL E. ; (17) 
"j 


and the scaled hydrogenic result for the contribution of this transition to 
D(u) is 
ro oxkT Pim 2 


Do. =—_= — —— a. foru>u, - Au. (18) 
ff Ry no. : 
G» YB j j 


i 


D(u) 


0 for u<u, - Au, ; 
J J 


where P; = fractional population of state i, 


m, = number of equivalent j-type electrons in state i, 


j 
a, = principal quantum number of electron j, 
u. = I./kT, 

oe 


Ry = Rydberg energy = 13.60 ev, 
@» = mean squared ionic charge in the plasma, 


Au, = reduction in 2 a owing to pressure ionization, 


and the free-free part of D(u) is 


2 
C sry? 1 764e 2 30 
Di? “a (ZH) ahwate ia, | (19) 


which is the result of taking the ionic charge as the effective charge for 


every free-free transition. (The bracketed quantity is 4. 76 x 10° cm“ /g. ) 
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We can determine the part of Die arising from free-free transition in the 
field of ions in state k: 
Peek te Cy 


where fe is just the net ionic charge of state k. Now the ratio P./P, is given 


by the Boltzmann-Saha equation: 


(20) 


e? j (21) 


where the g's are statistical weights. Hence, 


D. m, u . u,-Au, 

es ee Bi a | 
D = Ry n. 2 = (22) 
ff, k y 3, f. 8, 


So far, we have made no restrictions on the nature of the ionic states 
involved. We digress here and apply Eq. (22) to the determination of the 
absorption coefficient at low frequencies, so that large values of a are the 
only ones that contribute. For these states, since they are nearly hydro- 


. 2 
= = = 
genic, a = fe mi 1, and 8; /8, on Then 
2 
D, 2f-Ry u,-Au. 
a e? J for u>u.- Au. , (23) 
ff, k aT J J 


= 


and the part of D(u) due to all high bound states j on a core k, which we 


call D is given by 


bf, k’ 


bf,k j (24) 


where the summation covers only the j's for which u > ae - oy, > 0, since 
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the higher bound states have been pressure-ionized. Next, noting that 


du = 2i-Ry 
at npkT 


we replace the sum by an integral over us, recalling that in our approxi- 
mate pressure-ionization theory, au is independent of j for fixedk. We 


obtain 


D 
a a aoe (25) 
ff,k 


which implies for the whole continuous absorption coefficient that 


u 


D(u) = Die (26) 


out to some value of u at which the approximations based on large n start 


(9) (15) 


to break down. This result is given by Unsdld’”’ and by Raizer and is 


(4) 


also to be found in Eddington's book; we derive it here to show its 
invariance to the altering of the potential well implied by pressure ionization. 
Actually, this frequency dependence holds classically for any potential which 
deflects impinging electrons abruptly enough, since classically an electron 
radiates only while accelerated, so an instantaneous acceleration produces 
a uniform frequency spectrum of radiated power. Identifying this with the 
spontaneous emission coefficient and using the detailed balance relation 
Eq. (3), we obtain K~ Mg: 

Equation (26) is a surprisingly accurate representation of the mono- 
chromatic continuous absorption coefficients of air at temperatures of 


2 to 20 ev calculated by Armstrong, Holland, and Meyerott, its) 


generally 
deviating strongly only at values of u which are too large to affect seriously 
the Rosseland mean opacity. If Eq. (26) is inserted into the expression for 


the Rosseland mean, there results simply 


K = 1.15 D, 


R (27) 


f ’ 
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with Die given by Eq. (19). This provides a quick estimate of the Rosseland 


mean (neglecting scattering and line absorption), once G) is determined 


(15) 


for givenI and T. Raizer couples Eq. (27) with an approximate 


prescription for CG. but when a detailed occupation-number calculation 


is available, one may use the accurate Gry, When this is done using the 


(16) 


work of Armstrong, it is found that Eq. (27) reproduces his Rosseland 


mean opacities to ~30% at kT = 5 and 10 ev for all of his density range 


(o10-° to i078 afer): at kT = 2 and 20 ev the agreement is good at low 


densities but deteriorates toward higher densities, with Eq. (27) predicting 
opacities which are too large, but not by an order of magnitude. At very 


high temperatures, however, where u,, = 2*Ry/kT << 1, Eq. (22) shows 


K 
that bound-free absorption is negligible, and (still neglecting scattering) 


OE (28) 
KR ~ S(o) 196.5’ 


R by a factor of 226. As u, increases, 


Kp!/Deg increases--mostly because of K-shell absorption at first--and 


so that here Eq. (27) overestimates x 


reaches a value of 1/10 at u, = 1.8, so that Eq. (27) should be within an 


order of magnitude of the serach continuous opacity (without scattering) 
whenever U2 2. 

Returning now to Eq. (18), we may introduce an index, J, which labels 
a specific photoelectric transition, i.e., which specifies the initial state i 
and the initial quantum numbers m and e of the ejected electron. The 
total continuous D(u) is then obtained by summing over all possible photo- 


electric and free-free transitions: 


p,m 
Diu) = Dy [I+ RED EE ah ta - ys ay]. (29) 
(*> rT J 
where the unit step function H(x) = i" poe 3H takes care of the condition 
u>uy- Au,;, and the occupation-number calculation has set Py; = 0 whenever 


the initial state is pressure-ionized. Each term in the sum gives rise to 


an "edge, '' or jump, in the absorption coefficient. 
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The analogous expression of Armstrong, Holland, and Meyerott!!3) 
differs from Eq. (29) in two respects: each term carries a Gaunt factor, 


g(n,t .), and each term in which m, > 1 is replaced, in general, by several 


J 
terms weighted with fractional parentage coefficients which add up to m,. and 


depend on the quantum numbers Land S. The effect of the first siedeace 
is to redistribute strength among the edges corresponding to the same ae 
and different 4, leaving the total strength almost unchanged (since the 
weighted average over 4 of g(n£) is close to unity); and the second difference 
increases the number of edges, again without materially affecting the total 
strength. Test calculations of the continucus absorption coefficient of 
oxygen, using unit Gaunt factors and a level scheme in which each configu- 
ration is lumped at its center of gravity, have exhibited only minor differ - 
ences (~20%) in monochromatic absorption coefficients and negligible 
differences in mean opacities when compared with Armstrong's results. 

We suspect this is partly because the frequencies contributing heavily tothe 
Rosseland mean are around u ~ 5 and the transitions with m, > 1 tend to 


J 


have edges u, 2 @ (since the electron is being ejected from a substantially 


J 
populated level); and usually @ > 5. In our calculations we have therefore 
used configurations as our energy levels and neglected the Gaunt factors. 
The scattering contribution may be added at this point. We have 
taken it as the classical Thomson approximation to incoherent Compten 
scattering by free electrons. Since scattering by bound electrons occurs 
and is of the same order of magnitude, and since scattering dominates only 
when most of the electrons are free anyway, we have assumed that the 


Thomson cross section applies to all electrons, bound and free. The 


resulting contribution to D(u) is 


3 - 
u Z 8r 0 
Da = “4 
l-e 1 
: : Sid 3 2 
where r, = classical electron radius; the bracketed quantity is 0.398 cm /g. 


0 


V. EFFECT OF LINES ON LOCALLY AVERAGED 
ABSORPTION COEFFICIENTS 


In most cases, when line absorption has an important effect on 
frequency-averaged absorption coefficients such as the Rosseland and Planck 
mean opacities, the effect is due to a large number of lines scattered through 
the spectrum. Under such circumstances, the concept of a locally averaged 
absorption coefficient is pertinent. Bythis is meant an average with respect 
to frequency, taken over an interval large enough to contain several lines 
but small enough that the continuous absorption coefficient, Planck function, 
average line strength per unit frequency interval, etc., do not change sub- 
stantially within the interval. We will consider two such averages, the 


"local Rosseland mean" and "local Planck mean," defined respectively by 


Au 
Kplu) du' 
Au e(u') 
and 
it K(u')du' 
Au 
K (u) _ ’ 
e Ai du' 
Au 


where the integration range, Au, of the frequency variable u' is centered 
on u and satisfies the above conditions, and two models of the line absorp- 
tion coefficient. The statistical model, in which the lines are randomly 


located, has been worked out by Maver Elsasser's model, ne 


with 
uniformly spaced lines, is another case offering some analytic tractability. 
In this section we will compare the results of using these two models to 


compute the local mean absorption coefficients. 
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We consider a frequency interval in which the absorption coefficient 
consists of a continuous part, which we take to be constant and normalize 
to unity, and a line part denoted by r(u); we use the reduced photon frequency 
variable u = hv/kT. Let there be N identical lines in an interval Nd, with 
the center of the line at ui let a be the contribution of the ae line to r. 
Since all the lines are alike except for location, a depends only on u - uw 
Define the strength, S, of the ae line asf ry, du. The average line 
strength per unit frequency interval (in the present dimensionless units) is 
R=r-= S/d; hence, the local Planck mean is just 1 + R, which is independent 
of line shape. Schwarz’ inequality demonstrates that 1 + R is also an upper 
bound to the local Rosseland mean; it is the result of complete blending of 
the lines. A lower bound is, of course, unity, corresponding to infinitely 
sharp lines. 

To calculate the local Rosseland mean, it is necessary to specify the 
shape of the lines and some properties of their relative positions. Here 
we will consider lines with Lorentz (dispersion) profiles such as result from 


natural and collision broadening. Then 


w 


AWM 


a ; (u 7 +w 
“n 
where w is the half-intensity half-width in units of kT. In the statistical 
model it is assumed that the N values u(n = 1, 2, --: N) are randomly 
located in the frequency interval Nd, so that the joint probability that the 


center of the first line is in du, atu,, the second in du. atu etc., is just 


1 2 2’ 


N 
unayyN TT au, 
n 
n=1 
and the value of r(u) with this arrangement is ae a The average of 


1/(1 + r) over all arrangements is 


in ore co -rt 
« l > =¢nay NY rf P f es es jidu =a | dt a Te n du ; 
l+r u toon { n 
l+> or ie e ms 
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n 


29 


-(l+r)t 


where the fact that 1/(1 + r) 5° e dt has been used instead of 


Mayer's approximation 


4 -b.r 
i 


7 ™ D>. a 


i=] 


-rt 
Now all of the integrals [ e ™ du. are identical, so 


co -rt N foe) -rt 
GH «Joes be [eT a) -[Peet(geso-e bay) 


since the integral over u, covers the interval Nd. As N becomes large, 


1 
with d (the average line spacing) fixed, this integral becomes independent 


of N and u, and we get 


< 


co 
0 


-rit 
where G(t) - (af? a -e er : 


evaluated as 


For Lorentz profiles, G(t) can be 


G(t)= Rt(==); p=="; 


2 
r= 2 ity) 


iD 
~~ 
—_ 


) dy 


me [1() +1,(3)] 


1-5, x<< 1 


2 
Nix 


; x>> 1, 


and asa + r) depends only on the two parameters Rand f. Since in this 
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case the averaging over u has in effect already been carried out, we can 
identify as + r)» with the reciprocal of the local Rosseland mean, and 
denote it by T,(R, 8). 

is (3) 


In Mayer's report a further result obtained by Jacobsohn is presented: 
if the line strengths are not identical but have an exponential probability 
distribution with a mean value S, the resulting TAR, B) differs only in that 

the function F(x) must be replaced by 1/N1 +x, which is numerically quite 
similar. 


(17, 18) 


We turn now to Elsasser's model, in which the line centers are 


equally spaced, so that 


r(u) = eh — : 
Me ees (u - Hye + 7 


is periodic inu. As Elsasser showed, this can be summed into a very 
simple form: 


sinh 


Be cosh B - cos 0 


, 


where, a: before, R = S/d; B = 2mw/d; and 6 = 2mu/d. To average over u, 


we need only average over 6 from O tom. Thus, 
(ry eR. 
ea ae coth B 
r R : 


ane ee 


J R*4142Rcoth@ 


Comparing Tu with T_ for the same R and 8B, we find some interesting 


E 
results. For example, in the limit of small B, 


Bao PER ch be 
Ty =) /® +—BR+ ; 


fet n JOE AEE ac 
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and at large 6, 


iR 


2 
1 R \*1 
Tu cralit(sFe) zs? 


2 
l R -28 se] 
T, = peel +(e) ee + , 


so that the leading terms at extreme values of B are identical, but for . 


I 


intermediate B's the statistical array gives larger T, i.e., a smaller local 
Rosseland mean. Physically, this stems from the occurrence of occasional 
large ''windows" in the statistical array. 

Of perhaps more relevance to the present study is the fact that in 
either model the parameter 6 appears instead of w andd separately. This 
means that changes in the number of lines per unit frequency interval affect 
the opacity in much the same way as changes in line width, if the strength 
per unit frequency interval, R, is fixed. 

Cases sometimes occur in which the locally averaged absorption 
coefficient as used above is not meaningful. An example is iron ata 


(3) 


temperature of 1000 ev, which was treated by Mayer. A large part of 
the line effect comes from a fairly closely spaced group of strong 1s ~ 2p 
lines with centers ina fairly small frequency interval, and wings which 
extend far outside the interval. In the wing region, containing no line 
centers, the foregoing recipes would give no line contribution. Mayer used 
the statistical method in a form which correctly takes the wing effect into 
account; at distances from the group of centers large compared with the 
spread of the group, the wing absorption coefficient is independent of 
everything except the product of strength and width summed over all the 
lines in the group, i.e., the total strength times a strength-weighted average 
width. The effect of this group on the opacity is then insensitive to number 
of lines, but still depends on width, in contrast tothelocal averages, which 
depend on both in the same way. Generalizing, we may say that the change 
in the opacity produced by splitting each line in the spectrum intoN compo- 
nents shouldbe at most equalto the change produced by multiplying all 
widths by a factor N. 


VI. LINE BROADENING 


The theory of the broadening of lines in plasmas has developed consi- 
derably in recent years, partly owing to its importance in optical measure- 
ments of the electron densities of laboratory plasmas. A comprehensive 


(19) 


review by Baranger will shortly be available, and we extract from it 
the bulk of this section. 

Two sources of line width which persist to vanishing densities are 
Doppler broadening and natural broadening. However, both these effects 
are quite small in our range of interest. The Doppler width is ~NkT/AMc 
times the line frequency, and even in hydrogen at a temperature of 100 ev, 
this ratio is only ~3 x 10. The natural width is roughly (1/137)(Z/ 137)° 
times the spacing of lines in a series, and is negligible except possibly at 
very large Z. We have neglected both these broadening mechanisms. 

Pressure broadening is the name given to the line width produced by 
interactions between the radiator and other particles. In plasmas which 
are even slightly ionized, neutral perturbers may be neglected and the 
relevant interactions are Coulombic. Usually, most of the broadening is 
due to perturbers which are distant enough that the radiator feels only a 
homogeneous electric field which changes in time as the perturber passes 
by, and in all practical cases it is permissible to treat the perturber's 
motion classically. The spectrum radiated by the total system is then 
calculated quantum-mechanically. The theory simplifies considerably in 
the limiting cases of slow or fast perturber motion, which give rise to the 
"quasi-static approximation" and the "impact approximation, '' respectively. 

It turns out that the electron perturbers always obey the impact 
approximation in a plasma at nondegenerate densities (I >> 1). This 
approximation implies that the typical collision is weak and may be calcu- 


lated using perturbation theory; occasional strong collisions may also 
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occur. The result is that an isolated line acquires a Lorentz shape with a 
width proportional to the electron density, and a shift which is usually of 
the order of the width. When a closely spaced group of states gives rise 
to a group of overlapping lines, there are some asymmetric terms in the 
spectrum besides the Lorentz profiles. Most of the electron width is 
usually due to collision-induced transitions involving the upper state of the 
line, which give it a finite lifetime. 

Ion perturbers, on the other hand, can often be treated in the quasi- 
static approximation; this is true especially for the ions that give rise to 
the wings of the line profile. The effect of the ions is then a static Stark 
effect (which shifts levels adiabatically but induces no transitions), averaged 
over the probability distribution of the electric field produced at the radiator 
by the ions. The resulting profile depends on this probability distribution 
and on the Stark pattern of the line. 

The combined effect of ions and electrons is usually calculated by 
choosing an ionic field, finding the corresponding Stark pattern, broadening 
it by electron impacts, and then averaging over field strengths. Finally, 
possible corrections for ion motion are evaluated. This program has been 
carried out for a large number of lines of hydrogen and helium, and the 
resulting profiles are in excellent agreement with arc and shock-tube 
experiments, in which temperatures of 1 to 4 ev are produced. 

In our application we are mostly interested in the line absorption 
coefficient outside the line cores, in the far wings. Asymptotic wing profiles 
obtained using the full ion-electron broadening theory show that here the 
ion and electron contributions are additive, and that the ion part falls off 
more rapidly with distance from the line center. This is especially true 
when the radiator is itself anion, since the Coulomb repulsion between the 
ions inhibits the close encounters which prcduce the line wing. The electron- 
width calculation also is different in this case, since the classical perturber 
trajectories are hyperbolae instead of the straight lines which electrons 
follow near neutral atoms; somewhat surprisingly, the Coulomb attraction 


decreases the electron width slightly. 
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The evaluation of the electron width of lines emitted (or absorbed) by 
ions, and the accompanying asymmetric contributions, are given in detail 


by Baranger and Stewart. (20, 21) 


The width is proportional to the total 
inelastic cross section of the ion in the upper state of the line, and this in 
turn is given by perturbation theory in terms of the same dipole matrix 
elements as determine the line strengths, and of a Gaunt factor for 
bremsstrahlung which takes care of quantum corrections to the perturber's 
motion. When this is set equal to unity, a sum rule permits the evaluation 
of the cross section in terms of the mean square radius of the optical 
electron's wave function. In our problem, most of the lines will be high 
series members with upper states which are approximately hydrogenic, and 
the final result is thus 

2.0 


. n 
as aa a ce 


N 


where w = half-intensity half-width in units of kT, 

n = principal quantum number of opticalelectron inupper state, 

z = core charge of ion = net ionic charge plus 1. 
The total line absorption coefficient is composed of a Lorentz profile for 
each line in the spectrum, some asymmetric terms, and some ion-broadening 
contributions. The relative importance of these can be assessed by consi- 
dering a single series of lines having a common lower state and upper states 
which differ only inn. The asymmetric terms will tend to cancel each 
other at large n, where the line series becomes nearly periodic; they are 
relatively most important midway between line centers. We have evaluated 
them for some representative cases at low n, and find that even there they 
contribute only ~20% as much as the Lorentz part midway between centers. 
Likewise, the ion contribution can be overestimated by assuming linear 
Stark effect, and turns out to be negligible except very near the line cores, 


(22) 


where the absorption coefficient is high already. Griem has shown that 
for the higher series members, electron broadening dominates even in 


neutral hydrogen, which should be the case most favorable to ion broadening. 
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We have therefore kept only the Lorentz profiles, with half-widths 


determined as above. These widths are about half as large as those given 
by Sternheimer's result, which was used by Moszkowski and Meyerott'2?) 


(24) 


and by Cox. Though our expression was derived for lines emitted by 
ions, it predicts widths for neutral helium lines which are in qualitative 
agreement with the accurate electron widths calculated by Griem, Baranger, 


(25) 


Kolb, and Oertel; our widths are larger by a factor of generally less 


than 2. 


VII. LINE SERIES ABSORPTION COEFFICIENTS 


In this section we indicate how the line contribution is incorporated 
into the total absorption coefficient, using the pressure-broadened Lorentz 
profiles. The most straightforward (and longest) procedure would be to 
simply construct a list of all possible lines with their frequencies, oscil- 
lator strengths, widths, and initial-state populations, and at each frequency 
to find the sum of the absorption coefficients of all the profiles. This pro- 
cedure is fairly redundant because the line strengths are closely related 
to those of the photoelectric continua; the high series members at fre- 
quencies just below an edge frequency have an average oscillator strength 
per unit frequency interval which is the same as that of the continuum just 
above the edge. (This results from the correspondence principle.) In the 
same hydrogenlike, unit-Gaunt-factor approximation which we have used 
for the continua, the ae frequency dependence of the continuum extends 
also to the average strength per frequency interval of the whole line series. 

Using these facts, we can group the lines into series, with one series 
associated with each photoelectric transition. (Actually, there are in 
general two, since the optical electron can either increase or decrease 
its orbital angular momentum, /, by one unit; the lines of the second have 
nearly the same frequencies and widths as those of the first but are far 
weaker, and we ignore them.) Along such a series, labeling the lines by 
the principal quantum number, n, of the optical electron in the upper state, 
the line width grows as x whereas the line spacing shrinks. Eventually, 
the width exceeds the spacing and the series has coalesced to form a con- 
tinuum. Typically, this occurs at a lower frequency than that of the per- 
turbed edge (u, - Au,) which comes from the pressure-ionization correction 


to the edge position, so when lines are used the latter is superfluous. 
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(Of course, pressure ionization must still be included in the occupation - 


number calculation. ) 
We have assigned the line frequencies uw by assuming hydrogenic 


upper states, i.e., 


2 
f l 
Gye (32) 
cs aaa ene 
‘ n° kT 
where fy = net ionic charge of initial state of transition J, 
uy = unperturbed edge frequency of transition J. 

This does not imply hydrogenic line frequencies, since us is determined 
by an effective charge Z_, which in general is larger than f. + 1. 


J J 
We now introduce a closed approximation for the absorption coefficient 


due to an entire series. Recalling Elsasser's expression for a uniform 


band of lines 


r(u) 2 sinh B 


© cosh B - cos 6 (3) 


with 6 fixed and 0 linear in u, we replace B and 0 by smooth functions of 


u such that 


Blu) = anw /d , (34) 


O(u_) = 2m, (35) 


and use for qs (the local line spacing) the large-n approximation 


2 
2(f. + 1) 
gee Ue oe RY 
oh we kT ° (36) 


Inserting our previous expression for w. (Eq. (31) and eliminating 


n by Eq. (32), we obtain 


5/2 
B_(u) -2 (8 Re ar eer Ct ee (37) 


39 


1/2 
8 ,(u) = 20 (2x) (f5 + 1) (uy - ae et (38) 


and the contribution to D(u) of the sth photoelectric transition and its asso- 


ciated line series is now 


p,m 

_ Fok? *J J 2 

a aan G2) Ry 2 ae me) 
Zz J 

where the frequency-dependent part, which replaces the step function H of 


Eq. (29), is 
1 (u > u;) 


R (uy) . sinh B ,(u) j 40) 


cosh Bj(u)- cos Oyu) (43) 

This expression, which might be called an "osculating Elsasser band," 
generates Lorentz profiles when 8B is very small, overlapping lines when 
8 is fairly small, and a continuum when 6 is larger than about 1. The 
Lorentz profiles occur at the frequencies us and have the correct widths 
and strengths. Comparison of expressions like R (u) with actual sums of 
Lorentz profiles has revealed no serious distortion introduced by this 
approximation; all that occurs is a slight asymmetry in each line. 

One further point must be noted. At some low frequency 
uru,- ((¢,+1)’Ry]/[n5 kT] (where n, is the initial principal quantum 
number of the optical electron), R (u) will generate a fictitious same- shell 
line of grossly large oscillator strength (~100 in typical cases). The 
actual oscillator strengths of same-shell lines, when estimated using 
hydrogenic dipole matrix elements and actual line frequencies, turn out 
to be too weak to affect the opacity in our temperature range appreciably, 
since such lines occur at low frequency and have small widths. Accordingly, 


we discard the fictitious line by introducing a cutoff: 
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2 < Sse eee ee a 
R j(u) 0 for us<vy < ras 
ny +5 


(41) 


The position of the cutoff is reasonably irrelevant since in practical cases 
R (u) is usually small enough at the cutoff so that D(u) is dominated by 
other transitions and the free-free contribution. 

When the line contribution is represented in the way we have just 
described, the separation of the absorption coefficient into "line" and 
"continuous" parts is rather artificial. The "continuous" opacity could be 
defined as the result of placing each edge at the position uy - Au, (this is 
the usual convention), or at the position where B ,(u) = 1, or could be defined 
by some other criterion. In our results the ''continuous'" opacities we 
quote are the result of setting all the line widths equal to zero, so that 
the edges are at their unperturbed (zero-density) pcsitions: u,, andall 


J 
lines, including the fully merged ones, are omitted. 


Vill. ENERGY LEVELS 


The energy levels needed for the various stages of ionization and 
ae i ‘ S00 : : 26 
excitation are obtained by isoelectronic interpolation. Following Eagserc: 
we express the energy, relative to the bare nucleus with charge Z, ofa 


term y, ina particular configuration r, for an ion with N bound electronsas 


ae Z) _ W(x) 2 + Wi(roy, N)Z + Wilry N) + oz) : (42) 
where 
N l 
W,(r,N) = - De =. (43) 
i=] 2n, 


Wile, y, N) can be calculated from first-order perturbation theory and 
W (ry N) must be evaluated by comparison with experiment. Layzer has 


calculated the w,'s including interaction within a complex for all terms of 


the configurations Is "2s 2p". with q <2 andt< 6. 


A term-averaged W,, neglecting interactions within a complex, can 


1 
be expressed as 


g(r) 
W (rN) = 1% q(t) x ' (44) 
qinr n 
a 
where 
N = Sy Va (t) > (45) 
@inr 
eS De laule) = 8. lela eN (46) 
. Binr 4 Ba 
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and 
n* [a6] (a > 8) 
a 
1 2 
o(a,B) = 45 ny [aa] (a =68). (47) 
0 (a <8) 


The symbol q@ stands for a pair of quantum numbers (n, 4) according to the 
scheme a@(1s) = 1, a@(2s) = 2, a(2p) = 3, etc.; qay(*) is the number of @ 
electrons in the configuration r; and the quantity [a@B] is the two-electron 


(27) (11) 


interaction integral. Varsavsky lists the o(a@, 6)'s through a(4p) =8. 
Let E(r, f, Z) be the energy, relative to the gruund state of the neutral 
atom, of the configuration r in the "ta stage of ionization of an element with 


atomic number Z. Then, 


E(r, f, Z) _ : 5 2 
Ry [Wo(r,N = 28) Wi(r9.N = ZZ 
+ (W (r,N Zt) 5 Wy(ro.N 2. Z=f)) °Z 
Eq(r gf Z) 
+ [W (rN = Z-f) - Wolry,N = Z-f)] + —TRy (48) 


where ro is the ground term of the ion and Eo (ro f, Z) is the energy of To 
relative to the ground state of the neutral atom. 

In these calculations we have allowed as energy levels for each value 
of N: (1) all configurations in the lowest complex for N electrons with no 
more than one 3d electron; and (2) all configurations formed by adding one 
electron, bound by more than 2 ev if it is ann 24 electron, to all configu- 
rations in the lowest complex for (N-1) electrons with no 3d electrons. 

The term-averaged W's for those configurations covered by Layzer's 
calculations were used. For configurations with one or more electrons 
in the n = 3 shell, Eqs. (44), (45), and (46) were used to calculate the W's 


(11) 


using the o(a@,6)'s of Varsavsky. For configurations with an electron 
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with n24, shell-averaged energy levels were used. The W's for these 
configurations were calculated using Eq. (44) with shell-averaged c's. (3) 
All the W's were calculated by fitting to observed energy levels along each 


: ; 2 
isoelectronic sequence obtained from the NBS tables edited by Moore. (28) 


The error introduced by using isoelectronic extrapolation to obtain 


the energy levels is no greater than that introduced by using term-averaged 


levels. 


IX, CALCULATIONS 


The calculation is set up so that one need only specify the atomic 
number Z; the atomic weight A; a set of (Z + 1) energies [ the Eo(To: f, Z) 
for f = 0 to Z of Section VIII]; and the relative number of atoms of Z, Co 
for each element in a given mixture. For a given set of I and kT values, 
the equation of state, monochromatic absorption coefficients, and the various 
opacities for the mixture are then calculated at each I, kT point. The time 
required to perform the calculation at one ©, kT point on an IBM-7090, 
including the preparation on tape for plotting with an SC-4020 of the mono- 
chromatic absorption coefficients, is about 10 sec. 

The energy and statistical weight for each level and a complete 
specification of all possible transitions allowed between the levels are 
calculated, ordered, and placed on tape by a code called AUGEAS. It is 
not necessary to prepare a tape unless one wants to add an element or 


change the w)'s or W,'s in AUGEAS. Given an energy-level tape, one need 


only specify Z and C, for each element in the mixture and the I, kT values 
desired in order to carry out the calculation. 

AUGEAS contains a table of all allowed configurations with no n> 4 
electrons for N< 14. For each configuration there is a set of qy(r)'s 
(a = 1 to 6); a statistical weight Q(r); the quantities w,(r), w(t), and w/(¥); 
and an indication of whether or not attaching an n > 4 electron gives an 


allowed configuration. The quantities w,(r) are 


wi(r) = 2 (W.(r) - Wo(r9)] ; 
wir) =2[Wi(r)- Wilt]. (49) 


walt) = 2[ Wolr)- Wolr)] 
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where 


E(Z, r) = [w,(r)Z" - w (r)Z + wo(r)] Ry + E, (Fo) ; (50) 


The allowed configurations with one n> 4 electron are found by adding 
an electron to the indicated configurations in the basic table. For each 
stage of ionization f, n equals 4 to n max where n max is the next integer 
less than (f+ 1)VRy/2. The statistical weight for each of these configurations 
is 


Q = 2n“Qir) , (51) 


where r is the configuration to which the electron is added. The energy is 


2 
(Z - oc) 
a Ae 
n 


E = E(r, Z) - (52) 
where o is the screening of the n electron by the configuration r. 

For each Z, the energy and statistical weight of each configuration 
for f= 0 to Z are calculated. All possible transitions are obtained by 
matching configurations which differ by one in some qd, oF differ by not having 
an n > 4 electron attached to one of the configurations. The quantum number 
of the jumping electron, n, and the number of equivalent electrons that can 
make the jump, m, are obtained for each transition. 

The main code, DIAPHANOUS, uses the energy-level tape written 
by AUGEAS. An outline of the DIAPHANOUS calculation follows. 
I. Calculate for each value of f: 


1/3 


A. t-= [1+ 3(£+ DL] 1 
a 3/2 
Be teet 2 are fz 
(kT) TY 
and 

2 
f 

t. + 


fo 2. 
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The quantity Z* is an input number on the first pass. It is normally 
equal to one. Fi is an input factor that allows one to vary the amount of 


pressure ionization. It is normally equal to one. 


f-] 
eae 
B AE, = 7 J kT ; 
v=0 
AE, = 0 
Cc. X= 23.179 — 
(kT) 
3 (£4 1)3 
D. By, = Fy X 4, 128 x 10 ———s> 
(kT) 


X, and B; are used later in VI and VII. Fw is an input factor that 


allows one to vary the amount of line broadening. It is normally equal to 


one, 
II, Calculate for each energy level i: 


A. Set qa; = 0 for all AE. < Jy kT (QE, is the lowest transition energy 
i 
out of the state i). 


B. P Meee) 
: ,74,e 


e.=f4nT- 

ae 

€ = maximum of alle; : 
GP yo. OB 


all i for 
e.ement Z 
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CUP. 
p20 if Ec 107!* 
Zz 
III. Calculate: 
A, Ze > pf 
1 
Sn A 
B. Z => Pf 
i 
ae 
C. Z* = 2 ‘ 
Z 


* 
IV. Repeat I, HU, and Ill one time with Z_ as calculated on the first pass. 
V. CalcuPate: 


A. The electron density in number per cubic centimeter, 


3/2 
Ni. 26-052 07) 
e T 


B. The density in grams per cubic centimeter, 


Me 2. CrAr 


p = 1.673 x ioe ae 
Z 
C. The pressure in bars, 
P = 1,602 x 107 /8 € fk’ pep 
-/ e 


Zz 


D. The internal energy in ergs per gram, 


11 
9.56 x 10 2 _ 
E = 2 |= (1 + Z)kT + : 
SpA > ( ) ~ Pi fF, AEs] 
Z 
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E. The velocity corresponding to E, in km/sec, 


51/2 


V= 1.414 107 E! 


VI. Calculate for each transition j: 
A. u, = AE/kT. 
J J 


, 


B. $.= —t Pp. wu, 
J 2; Z J 


P, =P, in the initial state for transition j. 


2/7 
f. 
C. ulsu,-| 2 
j j on 
It is assumed that all lines are merged between a and ue 


D. An effective K edge is introduced for all levels in the ions with 


f<Z - 4. 
l, uy = 10.204 27 /kT =zsul. 
" “K : K 
2 
2. OK = 2PKUK , 
Z-4 
3. pL = ) P;. » 
* ie all i : 
for each f 
PK = 0 for Z< 4. 
4. fk = integer nearest to Zi 
5 ni = 1 
VII. Calculate: 
8 2° 


A. D = 80990 XK oO ———— SL 
ff 2 
[T (kT) a CA,] 
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B Ao = eee 3 ee : 
13. 605 22 
0. 977) C52 
Ge Res Con 
dS. C2Az 
zZ 
D. KRaizer = 1.15 Der : 


VII" In order to calculate KR: define a set of points: 


uy, = 0.55 +0.05k k = 0 to 30 
= 2.05 + 0. 025(k-30) = 30 to 230 
= 7.05 + 0. 05(k-230) = 230 to 330 
= 12.05+0. 10(k-330) = 330 to 360 
A. Calculate for each U k = 0 to 360. 
Wo a 
[ D(u, ) =D, e + Aol E(u,) + L(u, )] + K, =a, : 
l-e 
where 
Jk 
E(u, ) = De %; 
j=l 
and j, is such thatu, <u, <u', ; 
as ee 
J : 5 
sinh Bila) 


Liyy= 2) Clo, x(a), Blue 


. j cosh B.(u,) - cos x.(u, ) 
I-57 ‘ P; k jk 
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where 
x, < (2n,+1)7 , 
J J 
Cio.,x.,8.)=0 if<B.>2n, 
* 3 8;) B 
j-l 
o < 0.05 $4 ; 
J v=] " 
- “ufo 
B,= Blt - uy)", 
J 
-1/2 
x, = Xe (u, - uy) ; 
J 
and 
9 2 - ISgkT 
ar kT : 
2-5, kT 
= 0 if kT <0, 


: : ae 
where f' is the integer nearest to Z* ande_ is introduced to compensate 


for the edges between u = 2/kT and Jeu 


qi. 
ue 
oe SS = A(u,) = 3.849 x 1077 - =u 37 
u k ik, > 
[l-e °] 
(u, - u,) S (uy) 
] 0 v+] v-l 
3. = (u,) ” : H(u,) + 3 H(u) 
vzl 
(uy -u ) 
k k-1 
+ 7} H(u, ) ’ 
where 
A(u,) 
H(u ) = eons, 
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B. Calculate: 


] 
2 (u 


l. Kp = 


360) 
If the last edge occurs before u = 15.05, the remainder of the integral 


for «., is calculated with a series of five-point Gaussian quadratures. 


R 


3 Pp oie es 


2. pfu) = 3.849 x 10° ; 
(1 -e°")? [> (u + 0.05) - = (u - 0. 05)]} 


for u = 0.6 + n/10 (n= 0 to 144). This is the quantity that appears in the 


output as the absorption coefficient. 


[ pD(u!)] 
3. If any edges occur beyond u = 15, p(u!) = —- 
J u, 
J 
for all a > 15. 
2. In order to obtain kp: 


A. Calculate for each f: 


2 

. (f+ 1) 

1. X_ = 13.605 KT 

2 

; ~ x,/n Xr 
© Pee oe a 

n 

for 


n = 2 to n,max, 


n,max = the next integer less than [ 6. 803(f + 1)7] te ‘ 


x(n gmaxt 1)” 


“ngmaxtl,f” © ais 


B. Calculate for each transition j: 


Ng, maxt 1 
J 


—w 


n, f, 
a j 


ia 
C; KF 1.54 x 10 D,,[e + Ay o;e 


j R ; 
(1 + )! 
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X. DISCUSSION OF RESULTS 


The numerical results of our calculations are presented in Appen- 
dices A and B,which are self-explanatory. The range covered is 
kT = 1.5 to 34.0 ev; T= 20 to 8.4x io. Within this range the actual 
values used were chosen so as to give roughly uniform spacing in log T 
and log p. The density range depends on the temperature, but is roughly 
i107! to 107" vera: The elements we have treated are Z = 1, 4, 6, 7, 
13, 14; because the code gets its energy levels isoelectronically, it is 
immediately applicable to the intervening elements, and indeed to heavier 
elements provided that the temperature and density are such that virtually 
no ions may be expected to have more than 14 bound electrons. Thus, for 
example, iron (Z = 26) satisfies this condition if kT £n I > 355 ev, which 
is its thirteenth ionization potential. Provision for any mixture of such 
elements is also in the code. 

Appendix A gives the thermodynamic properties and mean opacities; 
Appendix B gives the 'tmonochromatic" absorption coefficients. To keep 
the presentation of the latter within reason, we have divided the range of 
u into two parts for each value of Z, I, and kT. In the first part, linear 
in u from 0.6 to 15.0, we plot the absorption coefficient p(= px), harmoni- 
cally averaged over intervals Au = 1/10. The second part is logarithmic 
in u from 15 to 1000 and shows the continuous absorption coefficient only. 
The values of p are scaled to the minimum values occurring in the regions 
u = 0to 15 and u = 15 to 150, respectively, and values more than 1000 times 
these minima are not plotted. Thus, the values of uw which contribute 
appreciably to the Rosseland mean are those which are exhibited in most 
detail. On each page, the pair of graphs appearing for u< 15 were calcu- 
lated with and without lines, showing the line effect on the ''monochromatic'' 
absorption coefficients. All of Appendix B includes the scattering contribution 


as part of the absorption coefficient. 
55 


The mean opacities, with and without lines, of Appendix A show some 
regularities in behavior which are exemplified (for the Rosseland mean) 
by Fig. 1. Though Fig. 1 was calculated for CH,: the pure-carbon or 
nitrogen results are not qualitatively different. The factor by which lines 
increase the Rosseland mean opacity, Kp (or decrease Ap) is seen to be 
correlated with the presence of strong photoelectric edges atu ~ 5, the 
important part of the Rosseland weighting function. Examination of Appen- 
dix A shows that this factor can amount to nearly 10 (hydrogen at 2.25 ev, 
high densities) and typically has a maximum in density (for a fixed temper- 
ature). The maximum comes about through a competition between avail- 
able line strength and line width. 

The effect of lines on the Planck mean, though much more dramatic 
(a factor ~10° at low densities), should actually be regarded as a caution 


against the indiscriminate use of « in emissivity problems. We recall 


from Section I that the Planck ee describes the net energy-loss rate of 
an object which is optically thin at all frequencies, including those of the 
centers of strong and narrow lines at Taege values of u; these lines are the 
ones responsible for the large line factors in the Planck mean, and in 
actual problems they will often be strongly self-absorbed, invalidating 

the use of the Planck mean as defined. (The validity condition is more 
stringent than just requiring the object to be small compared with a Planck 
mean free path.) Typically, Appendix A indicates a steady increase in the 
Planck line factor as the density decreases. 

The opacities without lines which we give, as mentioned before, have 
all the absorption edges at their zero-density (unperturbed) positions, so 
our line factors effectively come partly from the lines whose upper states 
are engulfed by the continuum due to pressure ionization. We have not 
attempted to give this part of the line factor separately. Our Planck mean 
opacities do not include scattering, whereas the Rosseland and monochro- 
matic results do (this is why Kp<Kkp in some cases). 


The absolute accuracy of our results is not easy to estimate, since 
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we think it is limited in most cases by the scaled hydrogenic matrix elements 
used. The strengths of individual edges may be in error by factors of 2 
or more, but compensation tends to occur, and seldom is a single edge 
crucial in the opacity. We can make some statements about the sensi- 
tivity to line widths: we have made a number of spot checks in which all 
line widths were doubled or halved. The effect of either of these operations 
is to change the Rosseland mean by about 10% when the line factor is about 3, 
and from the discussion of Section V it follows that the same or less effect 
should occur if we split each line into two components. A more detailed 
level scheme (such as distinguishing the separate LS terms of the L-shell 
configurations) would have as one effect such splitting, generally into 
several components. 

At lower temperatures or higher densities than those used in our 
calculations, effects which we have not included become important. We 
may mention degeneracy corrections, energy transport by free electrons, 
molecules, negative ions, and Rayleigh scattering. On the other hand, our 
method should be good up to temperatures of several kilovolts, where (at 


(29) 


low densities) the equilibrium positron density makes itself felt. 
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Appendix A 


THERMODYNAMIC PROPERTIES AND MEAN OPACITIES 


The results of our calculations are summarized in the following tables, 
grouped by element according to increasing atomic number, for each ele- 
ment grouped according to increasing temperature, and for each temperature 
listed according to increasing IT. There are two entries for each value of 
T: the first entry lists the opacities with lines; the second lists them with- 
out lines. 

The units of the temperature as printed are electron volts. The 
quantity Tis listed in the first column. The quantity in the last column, 


EGAM, is 


EGAM = 1 + 


wl 


The internal energy per unit mass, E, is listed under the column labeled 
ENG. That part of the internal energy due to ionization and excitation is 
labeled EION. The Rosseland mean opacity is labeled KROS. The Planck 
mean opacity is labeled KPLK. The numbers are printed in floating-point 


+ 
format; e.g., x.xxxx tyy is read as x. xxxx X 10 ed 
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